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Abstract 



A new approach is applied to the ID Anderson model by making use of a 
two-dimensional Hamiltonian map. For a weak disorder this approach allows 
for a simple derivation of correct expressions for the localization length both 
at the center and at the edge of the energy band, where standard perturbation 
theory fails. Approximate analytical expressions for strong disorder are also 
obtained. 

PACS numbers: 05.45+b,71.23.An,02.50.Ey 



I. INTRODUCTION 



Recently, it was suggested to treat ID tight-binding models with diagonal disorder in 
terms of classical Hamiltonian maps 0. This approach has been successfully used in the 
description of delocalized states in the so-called dimer model ||, as well as for the Kronig- 
Penney model ||. In this paper we show that even for the standard ID Anderson model, 
this approach allows us to obtain new analytical results and reproduce known results in a 
much more transparent way, by making reference to the properties of the dynamics of the 
noisy Hamiltonian map into which the model is transformed. 

As was indicated in |TJ,|2[], the discrete stationary Schrodinger equation 

4>n+l + 4>n-l = (Cn + E )^Pn , (1) 

with e n standing for the diagonal potential and E for the energy of an eigenstate, can be 
written in the form of a two-dimensional Hamiltonian map 

X n+ i = X n COS fJL-(j) n + A n X n ) Sin fi 

p n+ i = x n sin (j, + (p n + A n x n ) cos p . (2) 

Here, the variables (p n , x n ) play the role of the momentum and position of a linear oscillator 
subjected to linear periodic delta- kicks with the period T — 1. The amplitude A n of the 
kicks depends on time according to the relation A n = —e n / ship. For the Anderson model 
the distribution P(e) of the disorder is given by P(e) = 1/W for |e| < W/2, with variance 
(e 2 ) = a 2 = W 2 /12. Between two successive kicks, the rotation in the phase space is given 
by the eigenstate energy, E = 2 cos p. In such a representation, amplitudes ip n of a specific 
eigenstate at site n correspond to positions of the oscillator at times t n = n and, therefore, 
the structure of eigenstates can be studied by investigating the time- dependence of the tra- 
jectories in the phase space (p n , x n ). In particular, localized states correspond to unbounded 
trajectories and, vice-versa, extended states are represented by bounded trajectories. 

It is convenient to pass to action-angle variables (r n , 9 n ) according to the standard trans- 
formation, x = r sin 9, p = r cos 9. The corresponding map, therefore, has the form 

'"'ra+l ^nPn 

sin 9 n+ i = D^ 1 (sin(6» n - p) - A n sin 9 n sin p) (3) 
cos 9 n+ i = D^ 1 (cos(#„ - p) + A n sin 9 n cos p) , 

where 



D n = yjl + A n sin(2^ n ) + A\ sin 2 9 n . (4) 
The localization length I is defined by the standard relation 



N-l 



1 

n=0 



l- 1 = lim (A V In ) = (In ) , (5) 



where the overbar stays for time average and the brackets for the average over different 
disorder realizations. The contributions to can be splitted in two terms 
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The second term on the r.h.s. is negligible because it is the average of a bounded quantity. 
It becomes important only when also the first term is small, i.e. at the band edge p ~ 0. 
Thus, apart from this limit, the localization length can be evaluated from the map using 
only the dependence of the radius r n on discrete time. The ratio r n+ i/r n is a function only 
of the angle 9 n and not of the radius r n , thus the computation of the localization length 
implies just the average over the invariant measure p{9), which is an advantage with respect 
to transfer matrix methods. Moreover, since r n+1 /r n is positive, there is no need to work 
with complex quantities. 

In a direct analytical evaluation of ([5]) one can, therefore, write 

T 1 = J P(e) ln{D{e, 9)) p{6) d9de , (7) 

where P(e) is the density of the (uncorrelated) distribution of e n , and p{6) stands for the 
invariant measure of the one-dimensional map for the phase 6, see (|j). We use here the fact 
that p{6) does not depend on the specific sequence e n , but can depend on the moments of 
P(e), particularly on its second moment cr 2 (see below). As one can see, the main problem 
is in the expression for p(6), which was not found explicitly even in the limit of a weak 
disorder, A n -»• |]||. 



II. WEAK DISORDER 

By weak disorder we mean that A n is small. This can be arranged even at the band- 
edge, where the denominator sin p of A n is also small; thus the disorder e n must go to zero 
faster than p (how much faster, it is determined by the properties of the Hamiltonian map). 
Retaining only terms up to 0(A n ) in the map (|3p for 6 n one gets 

#n+i = #n ~ P — A n sin 2 n + A n sin 3 9 n cos 6 n mod 2n , (8) 

which coincides with formula (62) in Ref. J|. 

The expression ([?]) for I" 1 can be written in the weak disorder limit explicitly, 

r 1 = / e 2 P(e) de j p{6) Q - - cos(2£) + \ cos(4fl)) dff , (9) 

^ o 

which is valid over all the spectrum except at the band edge, where the additional contri- 
bution in (H|) is present (see below). In fact, standard perturbation theory [|] corresponds 
to the assumption that p(Q) is constant. Thus, one easily obtains 



~ Ssin^ 96(1- f) ' V ' 

This expression was found to work quite well over all energies, but, surprisingly, numerical 
experiments showed a 4% deviation at the band center. One had to explain why standard 
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perturbation theory fails at the band center, while it is correct everywhere else. Non- 
standard perturbation theory methods were devised in Refs. where the correct value 
for l~ x at the band center was obtained and, moreover, a different scaling with disorder 
was discovered at the band-edge ||. However, these methods hide behind mathematical 
difficulties the physical origin of the discrepancy. We are here able, by looking at the 
properties of map @, to understand both the physical nature of the discrepancy at the 
band center and the different scaling at the band edge. Moreover, our derivation is much 
simpler mathematically and more straightforward (this can be already seen from the very 
simple derivation of the expression ( [TDD ) . 



A. The band center 

In order to derive analytically the correct expression for l~ x at the band center E — , 
one has to find the exact expression for the invariant probability measure p{0). This latter 
arises from map (§) specialized to the value /i = ir/2. For vanishing disorder the trajectory 
is a period four, specified by the initial angle 9q. For a weak disorder any orbit diffuses 
around the period four, with an additional drift in 9. Asymptotically, any initial condition 
gives rise to the same invariant distribution, which can now be expected to be different from 
constant. To find this distribution, we write the fourth iterate of the map (|8|) 

2 

6 n+ 4 = n - sin 2 9 n - cos 2 9 n - °— sin (40 n ) , (11) 

where ffl = e n + e n+ 2 and ffl = e n+ i + e n+ 3 are uncorrelated random variables with zero 
mean and variance 2a 2 . Here, we have neglected in Eq. (|ll|) mixed terms of the kind e n e m 
(m n) and approximated {(fffl) 2 ) and ((ffl) 2 ) by their common variance 2cr 2 , which is 
meaningful in a perturbative calculation at first order in e n . 

Thus, the invariant distribution can be determined analytically in the continuum limit 
where 9 n+i — 6 n is replaced with dO and the random variables with the Wiener 

variables dW\ 1 dW2 with properties 

(dWi) = 
(dWidWj) = 25 ij a 2 dt i,j = 1,2 

obtaining the Ito equation 

2 

d9 = -dWi sin 2 9 - dW 2 cos 2 9 - y sin (49) dt . (12) 
To this we can associate the Fokker-Planck equation || 

7^(0>*) = y§9 ( sin ( 4 W^ *)) + [(3 + cos ( 4e )) p ^ *)1 • ( 13 ) 

The stationary solution p{9) of Eq. (|13|) , satisfying the conditions of periodicity p(0) = p(27r) 
and normalization J 27r p{9) = 1, is 
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p(9) = (W ^-Lj y/3 + cos(4fl) J , (14) 

where K is the complete elliptic integral of the first kind. One should note that the expression 
for the invariant measure has never been derived before, and it could turn out to be useful 
for obtaining observables other than Note that solution (|T^) does not depend on the 
strength of the random process. 

Inserting formula ([14]) into (|) at p, = tt/2 we get 



This result perfectly agrees with the one obtained in Ref . |5||| , although it is here derived 
with a different approach, that involves much simpler calculations. 

Thouless' standard perturbation theory result would correspond to neglect the average 
of the cos(46 l ) term in (|15|), meaning that the stationary solution ([14]) is approximated with 
a flat distribution. This approximation works well for all energies E = 2 cos fi, with p, = an 
and a irrational, but doesn't work for the band-center. Moreover, if one would consider 
observables which contain higher harmonics than those present in the formula for / _1 ([|), 
one would get corrections to standard perturbation theory also for other rational values 
of a = p/q. In fact, numerical experiments show that the invariant measure for rationals 
is modulated with the main period T = n/q (p and q being prime to each other). The 
amplitude of the modulation decreases with q; thus the strongest modification is obtained 
for a = 1/2, which corresponds to the band-center. It is moreover clear from formula (^|), 
that the only energy value for which a contribution due to the modulations in the measure 
p{9) is present in the inverse localization length I' 1 is the band center a = 1/2. This is 
because in Eq. (|J) only second and fourth order harmonics must be averaged, and only for 
a = 1/2 the fourth harmonic occurs in p(6). This is of course true only in the small disorder 
limit. 



B. The band edge 

The neighbourhood of the band edge corresponds to fi ~ 0. If the second order noisy 
term A 2 n in the map (H) is replaced by its average, which is the same approximation we did 
in the previous Section, the map (|8|) reduces to 

e n 5 2 

9 n +i = 9 n - At + — sin 2 9 n H — - sin 3 9 n cos 9 n mod 2n , (16) 
/j, fi z 

where 5 2 is the variance of the noise e n . For vanishing disorder and fi — > the orbits are fixed 
points. Moving away from the band edge produces a quasi-periodic motion and switching 
on the disorder gives rise to diffusion. Following the procedure of the previous Section (but 
here we do not have to go to the four-step map), we obtain the corresponding Fokker-Planck 
equation 

BP r) ( r) 2 \ r) 2 r) 2 

w {9,t) = - ^ - -sin 3 0cos0P(Mj + (sin^P(M) • (17) 
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There are in this case two small quantities: the noise e n and the distance from the band 
edge A = 2 — 2 cos /x ~ p 2 . Below we consider the double limit A — > 0, 5 2 — > 0. One can see 
that, if we keep the ratio k = p, 3 /S 2 fixed, the time scale of the drift term in ( |TTD is unique 
and, moreover, it coincides with the diffusion time scale, being We can thus rescale 

time t = tfj, and obtain the following stationary Fokker-Planck equation, 



^{k- sin 3 9 cos 9 p{9)) +~ (sin A 9 p{9)) = 



(18) 



which depends only on k . Its solution, with the same normalization and periodicity condi- 
tions as above, is 



P(0) 



sin 2 9 



C+ dx 



2J 



f(x) sin 2 x 



where 



f(9) = exp ( 2k ( | cot 3 6 + cot 9 



1 



(19) 



(20) 



and C, J are integration constants. To make p normalizable, constant C must vanish and J 
is then fixed by the normalization condition, 



J 



-l 



Stt 



(21) 



As was mentioned in Section I, in the evaluation of the inverse localization length given by 
(g) we must now take into account both terms on the r.h.s.; thus, we come to the expression 



in 



sin 9 



n+l 



sin 9 r 



(22) 



In the limit of a weak disorder and for /x — > one gets 



7T 

1 = -fi (cot 9 n ) = -2/2 J cot 9p{9) d9 . 



(23) 



After some straightforward calculations, with /i = (kd 2 ) 1 ^ 3 , one obtains 

_! (5 2 ) 1/3 f™dxx 1 / 2 exp(-£-2k 2 / 3 s 



2 / °° rfx x-V2 exp (- f - 2k 2 / 3 x 



(24) 



which coincides with expression (36) in Ref. 0. The limits k — > and k oo are then 
easily rederived and coincide with those in Ref. 0. For instance, the k — > limit gives the 
scaling law 



r J =^!^(5 2 ) 1 /3 = o.289... (5 2 ) 1 / 3 . 



2r (D 



(25) 
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It is interesting to observe that a similar scaling law was also found for chaotic billiards 
(stadia and oval ones) looking at the behavior of the Lyapunov exponent in the integrable 
limit Jl0[. It is quite natural to associate the Lyapunov exponent with the inverse local- 
ization length, and the geometrical parameter that in billiard measure the distance from 
integrability, with the intensity of the disorder v^ 2 in the Anderson model. 

We have seen in this Section that the study of the dynamics of the noisy circle map (§) 
allows us to derive both the standard Thouless perturbation theory results for the behavior 
of the localization length in the weak disorder limit, and the non-standard corrections to such 
a theory both at the band-center and at the band edge. The advantage of our method is that 
the derivation of the final formula has a clear physical meaning; moreover, for the band- 
center case, the procedure is mathematically more straightforward than those previously 
used@,§. 



III. STRONG DISORDER 

As is known, the analytical expression for the localization length was found only in the 
limiting cases of a very weak or a very strong disorder. It is interesting that relation (^) 
allows us to derive an approximate expression which is also good for a quite large disorder. 
Indeed, if the energy is not close to the band edge and the disorder is not very large, one 
can expect a strong rotation of the phase 9 . Therefore, the invariant measure p(6) can be 
approximately taken as constant, p{9) = (2ir)~ 1 . For such a disorder, one can explicitly 
integrate Eq.(^), first over the phase 9, 



1 2 r 1 / A 2 \ 

— J ln(l + Asin(2fl) + A 2 cos 2 #)d# = -ln 1 + — ; A 2 = e 2 /sin 2 /i 



(26) 



and after, over the disorder e, 



rj = \j m .„ U + *) d£ = I ,„ U + *L) + _ , (27) 

2 J V 4sm 2 /i/ 2 V 4sm 2 /i/ ^ 

Direct numerical simulations show that this expression gives quite a good agreement with 
the data for the disorder W < 1 -r- 3 inside the energy range \E\ < 1.85 . 

For a much stronger disorder, one can use another approach. Note, that for the unstable 
region 

\E - e n \ > 2 (28) 

of the one-step Hamiltonian map (§) both eigenvalues A^ 1,2 - 1 are real, 

1 



Ai ll2) = 2 {(E - e n ) ± y/(E - e n f - 4j (29) 

with A^A^ 2 -* = 1. Therefore, for stronger disorder W ^> 1 , one can compute the inverse 
localization length directly via the largest value A + of these two eigenvalues, by neglecting 
the region \E — e\ < 2, 
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-1 



(In |A+|) = | In Q (x + v 7 ^!)) dx = F{z x ) + F{z 2 ) (30) 



Here, x = \E — e\ and Z\ = W/2 + E, z 2 = W/2 — E; the function F is defined by 

F(z) = — (z In (z + Vz 2 - 4) - Vz 2 - 4 - z In 2) (31) 

Note that from this expression one can easily obtain the known expression for the localization 
length in the limit of a very strong disorder, 

W 

C^hiy-1 (32) 



IV. CONCLUDING REMARKS 

We have shown in Section II how to derive exact expression for the localization length in 
the weak disorder limit, using the properties of a noisy circle map (H). In particular, at the 
center of the energy band, a small correction to the localization length obtained by standard 
perturbation theory is needed, which is due to the contribution of the forth harmonics in 
the expression for the invariant measure p(6) . It is interesting that for other "resonant" 
values of the energy, the (weak) modulation of p(8) has no influence on the localization 
length. However, for other quantities than the localization length, these corrections may be 



important. In this sense, the exact expression (14) for the invariant measure p(9) obtained 
in this paper for a weak disorder, may find important applications. 

We would like to point out that the calculation in Ref. |j] of the localization length and 
of the density of states is performed after shifting the energy E — > E + x slightly away from 
the "resonant" values, the size of the shift being proportional to the variance x ~ a 2 for all 
energies except the band-edge. It is easy to see that the approach we have used here, also 
allows us for the derivation of the localization length near the center of the band. Moreover, 
we can also understand why the shift has to be, as in Ref. ||, of the order of the variance of 
the disorder. In fact, the recipe is that one should have A n ^> x, because, if the shift from 
a value of the energy corresponding to a rational value of a is too large, the orbit becomes 
quasi-periodic. In this case the modulation of the invariant measure which results in the 
non-trivial contribution to the localization length, is lost. 

Finally, it is interesting to note that the ID map (§|) can be compared with the Arnold 
map (llj 

8 n +i = 6 n — a + bsin 6 n mod 2n . (33) 

If we approximate A 2 n with its average, it is then tempting to associate the parameter a 
to our parameter p , and the parameter b to (A^). Although the modulation of the circle 
map (jS|) is a different function, and the noise is added through the term containing A n , our 
results show that the structure of the Arnold tongues persists (Arnold tongues are regions 
of the parameter space {a, b} where the dynamics is locked on a periodic orbit of period 
q, the tongues become narrower and narrower as b is reduced). Indeed, inside the tongue 
of the Arnold map any orbit corresponds to a rational rotation number p/q; outside the 
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tongues the motion is quasi-periodic. Trajectories inside the tongues of our model do display 
periodic motion with an additional diffusion, the periodic motion being responsible for the 
modulation of the invariant measure p{6). Outside the tongue, the motion is quasi-periodic 
also in our model and the invariant measure is flat. 
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